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Abstract 

Using the light cone QCD Sum Rules Method, we study the rare B uv^ decay 
and find that the branching ratios are, B{Bs vv^) ~ 7.5 x 10~^, B{Bd vu^) ~ 
4.2 X 10~^. A comparision of our results, on branching ratio, with constituent quark 
and pole dominance model predictions are presented. 
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1 Introduction 



The Flavour Changing Neutral Current (FCNC) process is one of the most promissing field 
for testing the Standard Model (SM) predictions at loop level and for for establishing new 
physics beyond that (for a review see and references therein). The rare decays provide 
a direct and reliable tool for extracting an information about the fundamental parameters 
of the Standard Model (SM), such as, the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements Vtd, Vts, Vtd and Kfe 

Right after the experimental observation of the 6 — 57 and B Xg'-y Q processes, 
the interest is focused on the other possible rare i?-meson decays, that are expected to 
get observed at future 5-meson factories and fixed target machines. Besides measuring 
the CKM matrix elements, the role played by the rare B-meson decays could be very 
important for extracting more information about some hadronic parameters, such as, the 
leptonic decays fs^ and fs^- Pure leptonic decays of the form, Bg and Bg 

are not useful for this purpose, since their helicities are suppressed and they have branching 
ratios 5(5, ~ 1.8 x IQ-^ and 5(5, l+l') ^ 4.2 x IQ-^^ §. For Ba meson case 

the situations gets worse due to the smaller CKM angle. Although the process 5, t^t^ , 
whose branching ratio in the SM is 5(5, t^t^) = 8 x 10^^ 0, is free of this suppression, 
its observability expected to be compatible with the branching ratio of the 5, 
decay, only when its efficiency is larger (better) than 10~^. 

Larger branching ratio is expected when a photon is emitted in addition to the lepton 
pair, with no helicity suppression. For that reason, the investigation of the Bs{d) l^l'j 
becomes intersting. Note that in the SM, the decay 5, —>■ vv is forbidden by the helicity 
conservation. However, similar to the 5, t'^t' case, the photon radiation process 
5, uu'-f takes place, without any helicity suppression. This decay is investigated in the 
SM, using the constituent quark and pole models as the alternative approaches, for the 
determination of the leptonic decay constants and Jb^ • shown in that work 

that the diagrams with photon radiation from light quarks give the dominant contribution 
to the decay amplitude that is inversely proportional to the constituent light quark mass. 
But the "constituent quark mass" itself is poorely understood. Therefore, any prediction 
in the framework of the above mentioned approaches on the branching ratios is strongly 
model dependent. Note that, similar obstacle exists for the 5,(fj) l^l'j decays as well 



In this work we investigate the Bg — * z/z/7 process, practically in a model independent 
way, namely, within the framework of the light cone QCD sum rules method (more about 
the method and its applications can be found in a recent review 0). The paper is organized 
as follows: In sect. 2 we give the relevant effective Hamiltonian for the b svv decay. In 
sect. 3 we derive the sum rules for the transition formfactors. Sect. 4 is devoted to the 
numerical analysis for the formfactors, where we calculate the differential and total decay 
width for the B z/z/7 ^ind confront our results with those of [^. Our calculations show 
that the constituent quark model and sum rules predictions are equal for the constituent 
quark mass ms{md) = (250) MeV. 

2 Effective Hamiltonian 

We start by considering the quark level process b qvv (g = s^d). This process is described 
by the box and Z-mediated penguin diagrams. The effective Hamiltonian for this process 
was calculated in ^ to yield 

^e// = Cg7^(l - 75)6z/7^(l - (1) 



where, 



2^2: 



a; + 2 3(x-2), , , 
X — \ (a; — 1)^ 



with. 



X 



In our calculations we shall neglect the QCD corrections to the coefficient A, since they are 
negligible (see for example 0). 

At quark level the process -Bs(ci) — * vv^ is described by the same diagrams as ^ gz/z/, 
in which photon is emitted from any charged particle. Incidentally, we should note the 
following pecularities of this process: 

a) when photon is emitted from internal charged particles (W and top quark), the above 
mentioned process will be suppressed by a factor (see 0), in comparision to the process 
b (lyy-i so that one can neglect the contributions of such diagrams. 

b) The Wilson coefficient C is the same for the processes b gz/z/7 and b — > qvv as a 
consequence of the extention of the Low's low energy theorem (for more detail see ||11||). 



So we have two types of diagram that give contributions to the process b — > qi'P^, when 
photon is emitted from initial b and hght quark hues. The corresponding matrix element 
for the process Bs(d) ^^1 is given as, 

{l\neff\B) = CD^.il - 75)i^(7l?7^(l - l^)m . (2) 

The matrix element (7|g7^(l — 75)6|i?) can be written in terms of the two gauge invariant 
and independent structures, namely. 



(7(9)197/^(1 -75)&|5(P + ?)) = V'47ra 



ml 



(3) 



Here, and g^^ stand for the polarization vector and momentum of the photon, p + q is 
the momentum of the B meson, g{p'^) and /(p^) correspond to parity conserved and parity 
violated formf actors for the B vPj decay. The main problem then, is to calculate the 
formfactors g{p^) and f{p'^)- For this aim we will utilize the light cone QCD sum rules 
method. 

3 QCD Sum rules for the transition formfactors 

Derivation of the effective Hamiltonian eq.(l) is one of the basic steps in the analysis of 
the Bq — >■ i/z/7 decay. We need to carry the calculation at the hadronic level, in another 
words, we must calculate the transition formfactors within the framework of some reliable 
theoretical scheme. We shall use QCD sum rules, more precisely, the light cone QCD sum 
rules method, to achieve this aim. 

According to the QCD sum rules ideology, one starts with the calculation of the transi- 
tion amplitude for the Bg — > uP^ decay, by writing the representation of a suitable correlator 
function in terms of hadron and quark-gluon parameters. So, to start with, we consider the 
following correlator: 

n^(p, q)-ij d^x e'^%^{q)\T [q^,{l - l,)b{x)bi^,q] |0) . (4) 

The general Lorentz decomposition of the above correlator is, 

n/.(p, q) = \/47ra [e^yapelPaqii^i + i [e* (pg) - (e*p)gj Hs} , (5) 
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with, Hi and 112 corresponding to the parity conserving and parity violating components 
of the correlator, and qp are the four-vector polarization and momentum of the photon, 
respectively. 

The formidable task here, is to calculate Hi and 112. This problem can be solved in the 
deep Euclidean region, where and {p+qf ^-re negative and large. The correlator function 
(4) in the framework of the light cone sum rules method was calculated in this deep region 
in |T^ (see also [|I3]-|T5[). We have recalculated this correlator and our final answer is in 



confirmation with the results of [Ij] . Omitting the details of the calculation, which can be 
found in [|l^, and after performing the Borel transformation for the formfactors g and /, 
the QCD sum rules method gives us: 



1 du 



u 
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(7) 



Here $('u)and gi{u) are the leading twist-2, while g^^^ and (7*^2) are the twist-4 photon wave 
functions, x is the magnetic susceptibilty, / = j-fpmp, with fp = 200 MeV |]12|, u = 1 — u, 
Cq and Cb are the charges of the light and beauty quarks, is the leptonic decay constant, 
and 5 = {ml — p'^)/{so — p^). The terms without the photon wave functions correspond to 
the perturbative contributions, when photon is emitted from heavy and light quark lines in 
the loop diagrams. The asymptotic form of the wave function is well known [P^]-||19||: 

= Quu . 
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The twist-4 wave functions entering in eqs.(6) and (7) are given by 



g(^\u) = --u{2 + u) 



4 Numerical Analysis 

The main issue concerning eqs.(6) and (7), are the determination of the gip^) and 
We first give a hst of the parameters entering in eqs.(6) and (7): 



= -(0.24 GeVf [19], (gg), = 0.8(gg), H, 
fe = 0.14 GeV [P , sq = 35 GeV'^ , rm, = 4.7 GeV 
\VtbV*\ = 0.045 , \VtdVZ\ = 0.010 m 



9p = 5.5 1^ 



(8) 



The value of x in the presence of external field was determined in p3|, ^ 

X{^^^ = 1 GeV'^) = -4.4 GeV 



r-2 



If we include the anomalous dimension of the current qcTa/sq, that is equal to — ^ at jU = m^, 
we get, = ''^ft) = —3.4 GeV^ . Following |T^, we shall take gi{u) = 1, to the leading 
twist accuracy. The Borel parameter M ^ has been varied in the region from 8 GeV^ < < 
20 GeV^. We have found that, within the variation limits of in this region, the results 
change by less than 8%. The sum rules for g{p'^) and /(p^) are meaningfuU in the region 
— ~ {few GeV"^), which is smaller than the maximal available value = m^. For an 
extension of the results to whole region of p^, we use the extrapolation formula. The best 
agreement is achieved with the dipole formulas (for more detail, see []12| and p5[). 

hi „, 2x h2 



9{P 



(1-4) 



2 ' 



with, 



hi ~ 1.0 GeV , mi ~ 5.6 GeV , 
h2 - 0.8 GeV , ma ~ 6.5 GeV 



Using eq.(2) and eq.(3) for the total decay rate, we get 



where, 

»i 



Here x = 1 — is the normahzed photon energy. Let us compare our results with the 



m| Jo 



(ix (1 — x)^ X |/^ (x) + (7^ . 



rriB 



ones that are obtained within the framework of the constituent quark and pole dominance 
models (Note that eqs.(6), (7), (15) and (16) in 0, are all misprinted and all these 
equations must be multiplied by the factor 3). The correct results are as follows: 



dx ml {AMf 



X (1 — x) 



(Constituent Quark Model) (10) 



dV _ C ag f^,mB,mB{l - xyx 
di " 1287r2 (m|. - xmlf 



1 



3 / I ZIT^ I ' (Pole Dominance Model) (11) 
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where, 



/ (y) = -17y' + A2y^ - 2% - 6 (4 - y) (1 - yf In {I - y) . 



The coupling constant for BgB*'~f transition in the constituent quark model is given by p6| , 

9 = +^. (12) 
ruq 



This coupling constant in the light cone QCD sum rules was calculated in |15| to give: 



^--T^- (13) 

jBjB*mB 

Using the values of the input parameters and the lifetimes t{Bs) = 1.34 x 10^^^ s, t{Bj) = 
1.50 X 10~^^ s [^, we calculated the branching ratios of the decays, Bs — > z/i^7 and 
Bfi VI"-). The results are presented in Table 1. The results in the Table for the third and 
fourth columns are obtained using the values for the coupling constant g given by eqs.(12) 
and (13), respectively. Note that, for the constituent masses, we used md — 0.35 GeV and 
m<j ^ 0.51 GeV . We find out that, eqs.(lO) and (11) yield results that are numerically close 
to eq.(9), with the constituent quark masses niq ~ fq\/2. If we set fq ~ 200 MeV we get 
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~ 250 MeV. If we use this value of the constituent quark mass, the branching ratios 
for Bg lyp^ and i^i^7, increase by a factor of 4 and 2.5, respectively. 

Also, for a comparision, we have calculated the photon spectra using the constituent 
quark, pole dominance and QCD sum rules models, and found that the photon spectra for 
the constituent quark and pole dominance models are fully symmetrical. But, as a result 
of the balance between a typical highly asymmetric resonance-type behaviour given by the 
non-perturbative contributions and a perturbative photon emission, the sum rules model 
yields a slightly asymmetrical prediction. 

In conclusion, we calculate the branching ratios for the processes Bg — > I'i^'y and Bd 
in SM within the framework of the light QCD sum rules and obtained that B[Bg — >■ 
uiyy) ~ 7.5 x 10"^ and B{Bd uu-f) ~ 4.2 x 10"^ Within this range of branching ratios, 
it is possible to detect these processes in the future B factories and LHC. 
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Sum rules 


Constituent 
quark model 


Pole 
dominance 


Pole 
dominance 




7.50 X 10-^ 


1.93 X 10-\^f 


1.79 X 10-%^ f 


0.94 X 10-^(^)2 




0.42 X 10-s 


2.26 X 10-9(^)2 


2.10 X 10-9(^)2 


0.52 X 10-9(^)2 



Table 
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